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HEAVY  TRAFFIC  RESULTS  FOR  SINGLE  SERVER  QUEUES 
WITH  DEPENDENT  (EARMA)  SERVICE 
AND  INTERARRIVAL  TIMES 


by 

P.  A.  Jacobs 

Naval  Postgraduate  School 
Monterey,  California 


0 .   ABSTRACT 
Models  are  given  for  sequences  of  correlated  exponential 
interarrival  and  service  times  for  a  single  server  queue.   These 
multivariate  exponential  models  are  formed  as  probabilistic  linear 
combinations  of  sequences  of  independent  exponential  random  vari- 
ables and  are  easy  to  generate  on  a  computer.   Limiting  results 
for  customer  waiting  time  under  heavy  traffic  conditions  are 
obtained  for  these  queues.   Heavy  traffic  results  are  useful  for 
analyzing  the  effect  of  correlated  interarrival  and  service  times 
in  queues  on  such  quantities  as  queue  length  and  customer  waiting 
time.   They  can  also  be  used  to  check  simulation  results. 


'Support  from  NSF  grants  ENG-77-09020  and  ENG-79-01438  and 
ONR  contract  NR-42-284  is  gratefully  acknowledged. 


1 .   INTRODUCTION 

Much  of  the  work  in  queueing  theory  deals  with  models  in 
which  assumptions  of  independent  service  and  interarrival  times 
are  made.   One  way  to  examine  the  effect  of  correlated  interarrival 
and  service  times  on  measures  of  performance  in  queues  is  to 
consider  models  in  which  the  interarrival  and  service  times  are 
correlated  stationary  sequences  of  random  variables  having  given 
marginal  distributions  and  compare  results  for  these  models  with 
those  for  which  the  interarrival  and  service  times  are  independent 
and  identically  distributed  with  the  same  marginal  distributions. 

In  [3]  Jacobs  and  Lewis  introduced  a  scheme  for  generating 

sequences  of  dependent  exponential  random  variables.   The 

EARMA  (1,1)  (exponential  mixed  autoregressive  moving  average  with 

both  autoregression  and  moving  average  of  order  1)  model  is  defined 

as  follows.   Let   (E  }   be  a  sequence  of  independent  exponential 

random  variables  with  positive  finite  mean   v   .   Let   {J  }   and 

{K  }   be  independent  sequences  of  independent  {0,l}-random  variables 

with   P{J  =1}  =1-3   and   P{K  =1}  =  1-p   where   0  <  3  <  1 
n  n  —    — 

and   0  <_   p  <  1   are  given  constants.   We  define 

(1.1)  X   =  3  E   +JA   ,,       n  =  1 ,  2  ,  .  .  . 

n      n     n  n-1 

where 

(1.2)  An  =  pAn_x  +  KREn  ,       n  =  1,2 

We  assume  that   An   is  an  independent  exponential  random  variable 
with  mean   v   .   Under  this  assumption   (X  }   is  a  stationary 
sequence  of  dependent  exponential  random  variables  with  moan  \> 


The  dependence  is  not  Markovian  in  general  and  as  was  shown  in  [3] 

(1.3)         Corr(X1,X1+n)  =  p11"1  (1-3  )  [  3  (1-p )  +  (l-3)p]  . 

Note  that  if   6=1   or   3=0   and   p  =  0,  then   (X  }   is  a 
sequence  of  independent  exponential  random  variables. 

The  EARMA (1,1)  process  is  formed  as  a  probabilistic  linear 
combination  of  independent  random  variables.   This  technique  can 
be  extended  to  form  sequences  of  multivariate  exponential  random 
variables.   This  idea  was  briefly  explored  in  Lewis  and  Shedler  [7]. 
These  models  have  the  advantage  that  the  marginal  distribution  and 
correlation  structure  of  the  sequence  are  specified  separately. 
In  section  2  we  will  give  several  ways  in  which  to  oenerate 
sequences  of  dependent  multivariate  exponential  random  variables 
from  probabilistic  linear  combinations  of  independent  exponentials. 
These  stationary  multivariate  EARMA  sequences  of  dependent  exponentials 
can  be  used  as  interarrival  and  service  times  in  queues.   Results 
for  measures  of  performance  for  these  queues  can  then  be  compared 
to  queues  with  stationary  independent  sequences  of  interarrival 
and  service  times  having  the  same  marginal  (exponential)  distri- 
butions.  Since  the  correlation  structure  in  the  multivariate 
EARMA  sequences  is  specified  by  parameters  independent  of  the 
marginal  distribution,  any  difference  between  results  for  queues 
with  these  multivariate  interarrival  and  service  sequences  and 
the  corresponding  queues  with  independent  interarrival  and  service 

uences  having  the  same  marginal  distributions  will  be  due  to 
the  correlation  structure. 

Suppose  that  the  interarrival  and  service  times  of  a  single 

ver  queue  form  stationary  multivariate  EARMA  sequences 
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of  dependent  exponential  random  variables.   Results  of  Loynes  [9] 
give  conditions  for  the  existence  of  limiting  distributions  for 
the  waiting  time  of  a  customer  and  the  queue  length.   However,  it 
seems  to  be  difficult  to  obtain  exact  results  even  in  the  simplest 
queues  with  EARMA  interarrival  and  service  times  (cf.  Jacobs  [4]) . 

As  a  result  we  turn  our  attention  to  the  possibilities 
of  obtaining  limiting  results  of  the  heavy  traffic  type  for 
queues  with  stationary  multivariate  EARMA  sequences  of  dependent 
interarrival  and  service  times.   Results  of  this  kind  will  indi- 
cate how  dependence  in  queues  introduced  by  the  EARMA  scheme 
affects  quantities  such  as  customer  waiting  time  and  queue  length, 
at  least  under  heavy  traffic  conditions.   These  limiting  results 
can  also  be  used  to  check  simulations. 

In  this  paper  we  will  show  that  Kingman's  [5]  result 
for  the  equilibrium  waiting  time   W   in  heavy  traffic  holds  for 
queues  having  multivariate  EARMA  interarrival  and  service  times; 
that  is,  we  will  show  under  conditions  detailed  in  Section  3,  that 
if  the  traffic  intensity   Ap     of  the  queue  is  close  to  1,  then 
the  distribution  of   (1  -  Ay   )W   is  approximately  exponential. 
The  mean  of  this  exponential  distribution  will  in  general  be 
different  from  the  case  of  independent  interarrival  and  service 
times.   Positive  cross-correlation  between  the  interarrival  and 
service  times  will  tend  to  decrease  the  mean;  while  positive 
correlation  within  the  interarrival  and  service  time  sequences 
will  tend  to  increase  it. 


In  Section  3  we  will  state  the  heavy  traffic  result  more 
precisely  and  present  several  examples  of  its  application.   Its 
formulation  for  multivariate  EARMA  queues  is  very  dependent  on 
the  fact  that  the  marginal  distribution  and  correlation  structure 
are  specified  independently  of  one  another  in  EARMA  formulations. 
In  Section  4  we  will  give  the  proof  of  the  heavy  traffic  result. 


2.   QUEUES  WITH  EARMA  INTERARRIVAL  AND  SERVICE  TIMES 

We  will  consider  a  single  server  queue  at  which  customers 

arrive  and  are  served  in  the  order  of  their  arrival.   The  Oth 

customer  arrives  at   t  =  0   and  finds  the  server  free.   The 

arrival  time,  service  time,  and  waiting  time  (excluding  service) 

of  the  nth  customer  are  denoted  respectively  by   T  =  ^k_n  X,  , 

S  ,  and   W  .   In  this  section  we  will  present  some  examples  of 
n         n  ^  ^ 

the  use  of  the  EARMA  scheme  to  model  dependent  interarrival  and 
service  times  in  queues. 

(2.1)   EXAMPLE.   Let   {X  }   be  an  EARMA (1,1)  process  with 

positive  finite  mean   A     and  parameters   0  <  3-.  <  1   and 

0  <  p,  <  1;  let   {S  }   be  an  EARMA (1,1)  process  independent  of 

{X  }   with  positive  finite  mean   p     and  parameters   0  <_  8?  <_  1 

and   0  <_  p   <  1.   Note  that  if   8   =  1 ,  8   =  1 ,  or   3X  =  1, 

13 2  =  0,  p2  =  0,  or   3   =  0 ,  P;L  =  0 ,  32  =  1 ,  or  B]_    =    0,  p±    =    0 

3p  =  0,  p_  =  0,  then  the  queue  reduces  to  an  M/M/l  queue. 

Other  special  cases  include  queues  with  Poisson  arrivals  and 

EARMA (1,1)  service  times  independent  of  the  arrival  process, 

and  queues  with  EARMA (1,1)  interarrival  times  and  independent 

exponential  service  times. 

The  sequence   { (X  ,S  )}   is  stationary  and   {X  }   and 
^  n   n  n 

{S  }   are  independent.   However  the   X  's  (respectively   S  's) 
n  ^  n        ^        J         n 

are  positively  correlated  with 


(2.2)         Cov(Xn,Xn+k)    =    A    2p\    1(1-31)  [31d-P1)  +   d-31)p1] 


respectively 


(2.3)  Cov(sn'Sn+k)    =    M    2p2    2(1-32)  [32(1-P2)  +  d-32)p2]) 

from    (1.2). 

1  2 

Let      F         (respectively      F    )       be    the      n-algebra    generated 

bv      X,  ,...,X         (respectively      S,,...,S    )       and      G    ,,        (respectively 
J         1  m  1  m  -m+k  J 
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G         )' be    the    a-alqebra    generated   by      X    ,.  ,    X    ,,,-,,...     (respectively 
=m+k;  m+k         m+K  +  1  ^  ^ 


=m+k 

S  ,.,  S  Ll,  ,  ,  ,  ...)•   Let   L  (F1)   (respectively   L  (G1,, 
m+k    m+k  +  1  =m  J      =m+k 


denote  the  collection  of  real-valued  functions  having  finite 

second  moment  that  are  measurable  with  respect  to   F    (respectivel' 

=m  ^ 

2ro+k>'i=1'2-       If      fi    {    L2(Em'       and      9i    «    L2<§m+k»'    i   =    I'*, 
then    there    is    a      kn       independent   of      f .  ,    g.,    i   =    1,    2,    such    that 

for      k    _>    k„ 


(2.4)     |e[  (fxf2) (g1g2) ]  "  E[fif2]  E[gi92]l 

=     |E[f1g1]     E[f2g2l     -    Eff-J     E[f2]     EfgJ     E[g2]| 

=     |E  [f 1g1 ]     (E[f2g2]     -    E[f2]     E[g2]} 

+    E[f2J     E[g2]     {E[figi]     -    E[f1]     E[gi]}| 


(2.5)  <    E[|figi|]     5pk/2    E[f2]1/2    E[g22]1/2 

+    E[|f2|]    E[|g2l]     5p*/2    E[f2]1/2    E[g2jl/2 

<    5[pf2    +    Pk/2]    E-Kf^)2]1'2    EHg^)2]1/2 


where  (2.4)  follows  from  the  independent  of   (x  )   and   ^S  ) 

1  n  n 

and  (2.5)  follows  from  (6.5)  of  Jacobs  and  Lewis  [3].   Hence, 

if   f   (respectively   g)  has  finite  second  moment  and  is  measur- 

1     2 
able  with  respect  to  the  product   o-algebra   |   *  F    (respectively 

1       ? 
G   ,  x  G  ,i ) ,  then  for   k   sufficiently  large 
=m+k    =m+k 

(2.6)   |E[fg]  -  E[f]  E[g]|  <_   5E[f2]1/2  E[g2]1/2  [P*/2  +  P*/2]. 


Hence,  the  sequence   { (X  ,S  )}   is  asymptotically  uncorrelated 
in  the  sense  of  Rosenblatt  [10]  and  is  cjj-mixing  in  the  sense  of 
Billingsley  [1 ] . 

(2.7)  EX/AMPLE.   In  Lewis  and  Shedler  [7]  one  of  the  schemes 

that  is  proposed  for  cross-correlating  the  arrival  process  and 

the  service  times  is  the  following.   Assume   {X  }   and  {E    } 

^  n  n 

are  independent  sequences  of  independent  exponential  random 
variables  with  positive  finite  means   A     and   p     respectively 

The  nth  interarrival  time  is  assumed  to  be   X    and  hence  the 

n 

arrival  process  is  Poisson  with  rate  \.   Let  the  nth  service  time 

(2.8)  S   =  SE   +  J  (AU_1B  ) 

n      n     n       n 


where 


(2.9)  B      =    pB      ,     +    K    X 

v  n  n-1  n   n 


and      3/       p,     {J    }       and      {K    }       are    as    in    (1.1)    and    (1.2)    so    that 

{S    }       and      (X    }       are    now   cross-correlated   sequences.      We    will 

1  n  nJ  ^ 

assume  that   B„   has  an  exponential  distribution  with  mean   A 

Since   E    and    B    are  independent,  {S  }   is  a  stationary 

n         n  ^  n  2 

sequence  of  exponential  random  variables  with  mean   y   .   If 
3=1,  then  the  queue  reduces  to  an  M/M/l  queue.   The   S  's 
can  be  shown  to  be  positively  correlated  with 

(2.10)  Cov(Sn,  SR+k)  =  y~2(l-3)   P  ,    k  =  1,2,...  . 

Further,  the  nth  service  time   S    is  positively  correlated  with 

the  interarrival  times   X  ,  X   ,,...,  X,   with 

n    n-i         -i- 


-1  k 

(2.11)  Cov(S  ,  X„  .)  =  (Ay)   (1-3) (1-p) P 

n    n  — k. 


for   k  =  0 ,  .  .  .  ,n-l . 

Note  that   (X,S,...,X,,,S,,)   is  a  function  of 
n   n       n+k   n+k 

(B   t^X  ,  fc  ,  J  ,K  ,X   ,  ,  fc   ,  ,  J   ,-,,K  ,,,...,  X   ,,t   y    i  o        ,,K  ,1) 
n-1   n   n   n   n   n  +  1   n  +  1   n  +  1   n+1        n  +  k   n+k   n+k   n+k 

Hence,  since  the  sequence   {B  }   is  stationary  so  is  the 
sequence   {(X  ,S  )  .   Further,  by  arguments  similar  to  those 
used  in  obtaining  (6.5)  of  Jacobs  and  Lewis  [  3]  ,  if   f   and 
g   are  as  in  (2.6) ,  then  for   k   sufficiently  large 

(2.12)        |E[fg]  -  E[f]  E[g]|  <  5(>k/2  E[f2]1/2  E[g2]1/2 

so  that  the  sequence   {(X  ,S  ) }   is  A-mixinq. 

n   n 


There    are    many  ways    to   cross-correlate    the    interarrival 
and   service    times    using    the    EARMA   scheme   beyond    those    given 
above    and    in    Lewis    and   Shedler    [7] .      We   will    present   one   more, 
somewhat  more    general    than    that   of   Example     (2.7)     in   which    the 
service    times   were    coupled   back   over    the    previous    interarrival 
times  . 


(2.13)  EXAMPLE.  Let  {E  },  {E  },  and  {Cn}  be  independent 
sequences  of  independent  exponential  random  variables  having 
respective  positive  finite  means  X       ,    \i       ,    and  1.   Let 


(2.14)  Xn  =  6En  +  Jn(A-\) 


(2.15)  Sn  =  yEn    +  In(M_1An) 


with 


<2-16>  An  =  pAn-l  +  KnCn 

where   3,  p,  {J  }   and   {K  }   are  as  in  (1.1)  and  (1.2)  and  {1^} 
is  a  sequence  of  independent  binary  random  variables  independent 
of  everything  with   P{I  =1}  =  1  -  P{ln=0}  =  (1-y)   for  some 
fixed  constant   0  <_  y    <_  1 .   Again  we  assume  that   AQ   is  an 
independent  random  variable  having  an  exponential  distribution 
with  mean  1.   Then   {(En,  Er,  Ar  ,  In ,  Jn) 1   is  a  stationary 
sequence  of  random  variables  and,  hence,  the  same  is  true  of 

f<xn,sn)}. 


If   y  =  If   3=1   or   Y  =  1/   3=0,   p  =  0   or  y    =    0, 
3=1,   p  =  0,  then  the  queue  reduces  to  an  M/M/l  queue.   If 
3=1/   p  =  0   then  the  queue  reduces  to  a  special  case  of  Example 
(2.7).   Straightforward  manipulation  shows  that  in  general  the 
interarrival  times   X    are  positively  correlated  with 

(2.17)     Cov(Xn,  Xn+k)  =  A"2(l-3)2pk  ,        k=l,2,.... 


The  service  times   S    are  also  positively  correlated  with 

n  tr  j 


(2.18)     Cov(Sn,  SR+k)  =  y  2(1-Y)2pk,         k  =  1,2,... 


The  interarrival  times  and  service  times  are  positively  cross- 
correlated  with 


(2.19)     Cov(Xn,  Sn+k)  =  (1-8) (1-Y) (Ay)  1  pk ,  k  =  0,1, 


and 


(2.20)     Cov(Xn+k,  Sn)  =  (1-3)  (1-Y)  (Ay)  1  pk ,    k  =  0,1, 


Finally,  by  the  arguments  used  in  obtaining  (6.5)  in 
Jacobs  and  Lewis  [3]  ,  if   f   and   g   are  as  in  (2.6)  ,  then 
for   k   sufficiently  large 
(2.21)     |E[fg]  -  E[f]  E[g]|  <  5pk/2  E[f2]1/2  E[g2]1/2  . 


so  that  the  sequence   {  (X  ,S  )}   is  cb-mixina 

n   n 
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3.       A    HEAVY    TRAFFIC    RESULT    FOR   EARMA    QUEUES 

Let      {X    }        (respectively      {S    } )       be    the    sequence    of 

interarrival    times     (respectively    service    times)     for   a    single 

server   queue    with    first-in-first-out    service    discipline. 

Put      U      =    S      -    X        and   assume    that      {U    }       is    a    stationary 
n  n  n  n  J 

sequence   with      E [U    ]     <    0.       Let      W        be    the   waiting    time    of    the 
^  n  n  ^ 

nth    customer   and      W   =    sup     .   -.    ),     ,     U,        (with    the    convention 

^n  ^  0  /-k=l   k 

v0 
that   2,-1  =    0)  .       Loynes  [9]  has  shown  that  if  the  sequence 

{U  }   is  metrically  transitive,  then  the  distribution  of   W 
n  n 

tends  in  the  limit  to  the  distribution  function  of   W.   Since 

it  is  usually  difficult  to  obtain  the  distribution  of   W 

analytically,  particularly  when  the  service  and  interarrival 

times  are  dependent,  we  will  obtain  a  heavy  traffic  limiting 

result  for  it.   In  this  section  we  will  give  the  statement 

of  the  result  as  well  as  some  examples  of  its  use  leaving  the 

proof  to  the  next  section. 

For  each    a  >  0   let   Q    be  a  single  server  queue 

with  FIFO  queueing  discipline  having  stationary  sequences  of 

interarrival  times   {X  (a) }   with  positive  finite  mean   A (a) 

and  service  times   {S  (a)  }   with  positive  finite  mean   y  (a) 

n 

such  that,  if   U  (a)  =  S  (a)  -  X  (a),  then   {U  (a)}   is  a 
'       n        n        n  n 

stationary  sequence  of  random  variables.   Let 


(3.1)    E[U  (a)]  =  v(a)  <  0;  Var  [U  (a)]  =  s  (a)  ■ 
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(3.2)     a  (a)  =  Var  U±(a)    +2  I      Cov (l^  (a) ,U1+n (a)  )  ; 

n=l 


_2 
(3.3)      d(a)  =  |v(a)|/a(a)     and    n(a)  =  d(a) 


We  will  assume  the  following 


(3.4)   ASSUMPTION.   There  exists  a  distribution  function   G   so 

2 
that   lim      P{[U  (a)  -  v(a)]   <  x}  =  G(x)   for  all  continuity 

a  ->■  °°      n  —  J 

2       2 
points   x   of   G;  lim      v(a)  =  0;  and  0  <  lim    s  (a)  =  s  <  <*>. 


(3.5)  ASSUMPTION.   0<lim   n2(a)  =  a2  <  <~  . 

(3.6)  ASSUMPTION.   If  f  (respectively   g)  is  a  function  with 
finite  second  moment  measurable  with  respect  to  the  o-algebra 
generated  by   ^  (a)  ,  .  .  .  ,  U^  (a)   (respectively   Um+k(a),  um+k  +  1  (ex)  ,  .  .  .  )  , 
then 

(3.7)  |E[fg]  -  E[f]  E[g]|  <  <|>(k)  E[f2]1/2  E[g2]1/2 

for  some  decreasing  function  <p      bounded  above  by  1  not  depending 
on   ex   satisfying 


(3.8)  I       <j>(n)  < 

n=l 


Assumption  (3.4)  concerns  the  convergence  of  the  distribution  of  U2 (a) 

2  n 

[S  (.)-X  (<)]   and  its  first  two  moments  as   a  ■>  °°.   Assumption 


(3.6)  is  roughly  a  condition  that  the  sequences   (U  (a),  n=l,2,...} 

12 


be  <j)-mixing  uniformly  over  all   a.   This  assumption  guarantees 
that  the  series  in  (3.2)  converges  and 


i      n  2      2 

lim   -  E[(  I       (U,  (a)  -  v(a)  ))  ]  =  a  (a) 

n  -*■  °°        k=l 


Assumption  (3.5)  concerns  the  convergence  of  this  normalized 
variance  term  as   a  -*-   °°. 


Put 

[n(a)t] 

a  ■  k=i 


Z..(t)  =  d(a)  a(a)  X     I_    UR(a)  ,     t  >  0 


where   [x]   is  the  largest  integer  less  than  or  equal  to   x. 
Let  ID([0,«))   be  the  collection  of  real-valued,  right  continuous 
functions  on   [0,°°)   which  have  left  limits  everywhere  endowed 
with  the  metric  given  in  Lindvall  [8]. 


(3.9)   THEOREM.   If  (3.4)-(3.6)  are  satisfied,  then  the  process 

Za  =  ^Za^;t  -  °^   converges  weakly  in   B([0,°°))    to  a 
Brownian  motion  with  negative  unit  drift. 

The  proof  of  this  result  will  be  given  in  the  next  section 

Let 


W(a)  =  sup    I      U,,  (a) 


n  >_  0   k=l 
Then  it  follows  from  Theorem  (3.9)  that 
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lim      p{d(a)    w(a)     <    A  =      lim      P{sup     Z    (t)     <    x} 
a   -   oo   '     ^a(a)  "      J         a   +   -        t>  0 

=   1    -   exp(-2x} 

where    the    last   equality    follows    from   the    continuous   mapping 
theorem    (cf.    Billingsley     [1]).       Thus   we   have    the    following 
result. 

(3.10)       THEOREM.       If     (3.4) -(3.6)    hold,    then 


_  o 

lim      P{|v(a)|    W(u)     <    x)    =    1    -   exp{-2a      x) 


Hence,  if  the  traffic  intensity  Xy  of  a  queue  is  less 
than  but  close  to  1  and  (U  }  is  a  stationary  sequence  satisfy- 
ing (3.7)  and  (3.8),  then  the  distribution  of   W   is  approximately 

12       -1 
exponential  with  mean   -  y  a   E [U, ]     where 

oo 

a2  =  VartUjJ  +  2   J   Cov(ui'ui+n)  - 

n=l 


Thus  positive  cross  correlation  between  the  interarrival  and  service 
times  will  tend  to  decrease  the  average  mean  waiting  time;  while 
positive  correlation  within  the  interarrival  time  or  service 
time  sequence  will  tend  to  increase  the  mean  waiting  time. 
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(2.11)   EXAMPLE.   For  each   a  >  0,  let   Q    be  a  queue  as  in 
Example  (2.1)  with   {X  (a)}   (respectively   {S  (a)})   forming 
an  EARMA(1,1)  process  with  parameters   A  (a)  ,  3-,   and   p, 

(respectively   y  (a)  ,  3?   and   p~).   Assume   A (a)  <  y  (a)   for 

each  a,    but   lim      A (a)  =  lim      y (a)  =  m;  then  hypotheses 

a  -*00  ot-*   • 

(3.4)-(3.6)     are    satisfied   and   Theorem    (3.10)    holds   with 


a2    =    2m   2{1    +    (1-P;L)    X    (1-31)  [d-P1)31    +    (1-3]_)p1] 
+    (l-p2)_1     (l-32)  [(1-P2) 32    +    (1-32)p2]} 


Hence,  if  the  traffic  intensity   Ay     of  a  queue  with  independent 
EARMA(l,l)  services  and  interarrival  times  is  less  than  but 

close  to  1,  then  the  distribution  of  the  equilibrium  waiting 

1   2 
time   W   is  approximately  exponential  with  mean   ~  t  °   E  [U,  ] 

-1     -1         2 
where   E[U,  ]  =  y    -  A     and   a    is  given  by 

a2  =  A-2  {1  +  2(l-p1)"1  (1-31)  [(l-p1)31  +  (1-31)P2]) 

+  y"2{l  +  2(l-p2)"1  d-S2)  [(l-p2)32  +  (1-32)P2])  • 

Thus  the  dependence  in  this  queue  tends  to  increase  the  mean 
customer  waiting  time  at  least  in  close  to  heavy  traffic 
conditions . 


i  c, 


(3.12)   EXAMPLE.   For  each  a  >  0,  let   Q    be  a  queue  as  in 

Example  (2.7)  with  parameters  E[X  (a)]  =  A  (a)-1,  3,  and   p 

and   E[S  (a)]  =  y  (a)   .   Again  we  will  assume   A (a)  <  y(a) 

for  each  a   but   lim     A(a)  =  A  =  lim      y(a)  =  y  =  m.   Hence 

(J£-7"   OO  Q      ->.   OO 

(3.4)-(3.6)     are    satisfied   and   Theorem    (3.10)    holds   with 


o2    =    A"2    +    y"2  -  2(1-3)  (1  -  p)  (Ay)"1 

+    2p(l-p)~1     [(l-3)2y"2-    (1-3)  (1  -  p)  (Ap)-1] 
+    y"2[l    +    2(l-3)2    p(l-p)"1]  -  2  (Ay)"1     (1-3) 

=    2m"2  [1    +     (1-3)2    p(l-p)"1    -    (1-6)]     . 

2     2 
The  quantity      -  a .  ,   can  be  either  positive  or  negative  where 

7         —2—2—2  —1 

+  y    =  2m   .    It  is  negative  if   p  <  (2-3) 

ind 

and  nonnegative  otherwise.   Hence,  if   p   is  sufficiently  less 
than    (2-3)   ,  then  the  average  customer  waiting  time  in  the 
queue  will  tend  to  be  less  than  the  waiting  time  in  the  inde- 
pendent case  at  least  in  close  to  heavy  traffic  conditions. 
It  will  tend  to  be  greater  if   p   is  sufficiently  greater  than 
(2-3)"1. 


(3.13)   EXAMPLE.   For  each   a  >  0,  let   Q    be  a  queue  as  in 

a 

Example  (2.13)  with  parameters   A  (a),  y(a),  p,  3 , and   y. 
Assume   lim     ' (  )  =  A  =  lim      y(a)  =  y  =  m   as  before. 
Then  (3.4)-(3.6)  are  satisfied  and  Theorem  (3.10)  holds  with 
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o2    =    p2    +    X2-    2(Ay)  1     (1-3)  (1-Y) 

+  2p(l-p)"1(,r1  (1-y)  -A_1(l-3)}2 


=  2m  2[l-(l-3)  (1-y)  +  Pd-P)  1{(3-y)  }2] 


2     2 

The  difference   a   -  a.  ,   can  be  either  positive  or  negative; 

2       -2     -2      -2 

(a.     ,  =  A    +  u    =  2m   ) .   The  difference  will  be  neqative  if 
ind 

2 
p  <  (1-3)  (1-y)/ [(1-3)  (1-y)  +  (3-y)  ]   and  nonnegative  otherwise. 

Hence  under  close  to  heavy  traffic  conditions,  the  average 

customer  waiting  time  will  tend  to  be  less  than  the  independent 

2 
case  if   (1-3)  (1-y) /[ (1-3)  (1-y)  +  (3-y)  ]   is  sufficiently  greater 

than   p.   It  will  be  greater  than  the  independent  case  if  the 

same  quantity  is  sufficiently  less  than   p. 
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4 .   PROOF  OF  THE  MAIN  RESULT 

In  this  section  we  will  give  a  proof  of  Theorem  (3.9) 
We  will  assume  that  (3.4) -(3.6)  are  satisfied.   Let 


,,  ,  [n(a)t] 
Ya(t)  =  ^|     ^    [Uk(a)  -  v(a)] 


By  Theorem  31  of  Lindvall  [8]  and  the  continuous  mapping  theorem, 
Theorem  (3.9)  will  be  proved  if  for  each   a  >  0   we  can  show 
the  following  result. 


(4.1)   THEOREM.   The  process   {Y  (t);  0  <_  t    <_   a}   converges 
weakly  on  D([0,a])   to  Brownian  motion. 


(D([0,aj)   is  the  space  of  real-valued  right  continuous  functions 
on   [0,a]   with  left-hand  limits  endowed  with  the  Skorohod   J, 
topology. ) 

We  will  prove  Theorem  (4.1)  for   a  =  1   by  a  series  of 
lemmas.   The  proof  for  arbitrary   a  >  0   is  similar. 


(4.2)   LEMMA.   The  process   {Y  (t) ;  0  <  t  <  1}  has  asymptotically 

a        —   — 

independent  increments. 


Proof.   Let   0  <  sn  <  tn  <  sn  <  tn  <  *  *  •  <  s   <t   <  1.   Let 

1—1     2—2  r  —   r  — 

C.   be  the  event   {Y  (t.)  -  Y  (s.)  <  x.}  for  real   x..   Then 

l  a   l      a   l   —   i  l 

C.   is  an  event  in  the  o-algebra  generated  by 
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U[n(a)s.]+1'  •••  '  U[n(a)t.r   If   6   is  the  smallest  difference 


i  ~  ti-l'  then   tn<a)si]  +  1  -  [n(a)t._1]  >  [n(a)6].   By  (3.7) 


s  .  - 


r         r 
|P(  H   C .  )  -   n   P(C.  )  I  <  r<j>(  [n(a)  6]) 
i=l       i=l    1 


Since   6  >  0,  r<J)  (  [n  (a)  <S  ]  )   tends  to  zero  as   u  -*■   °°      and  the 
result  follows. 

In  the  following  results,  let   D. (a)  =  U. (a)  -  v (a) 
and    V  (a)  =  ).  ,  D.(a).   By  (3.4)-(3.6)  the  following  result 
holds . 


(4.3)  LEMMA 


oo 

a)  i  E[V,  (a)2]  <  2   J   |E[D, -(a)  D ,,.  (a)]|; 
K     K  _.=0      l      1  +  3 

b)  lim   n(a)"1  E[Vn(a)(a)2]  =  o2  • 

a   ■>  oo 

4        4      2 

c)  If   D,   is  bounded  by   C,  then   E[V  ,  .  (a)  ]  <_  K  C   n  (a) 

where   K    depends  on  $      alone. 


(4.4)  LEMMA.   The  sequence  of  random  variables 

-1  2 

{n(a)   (V  ,  v  (a))  ;  a  >  0}   is  uniformly  inteqrable 
n  (a)  — 


2 
Proof.   By  (3.4)  ,  {D-.  (a)  ;  a  >  0}   is  uniformly  integrable.   The 

remainder  of  the  proof  is  very  similar  to  that  of  (20.48)  on 

page  176  of  Billingsley  [1]  and  will  be  omitted. 
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(4.5)  LEMMA.   For  each  positive   £   there  exists  a   y,  y    >  a 
and  an  integer   nn   such  that   n(a)  >_  n    implies 


P{  max    |vi(a)  |  _>  y/n  (a)  }  <  -^ 

i  <  n  (a)   1  y 


Proof.   Put   C  (a)  =  J1  ,  I D . (a) | .   By  (3.4),  there  exists 

an  increasing  sequence  of  integers   m.   such  that   n (a)  >  m. 

2  2  -1 

implies   P{|D-.(a)|  >_  y    /n  (a)  /i  }  £  (y  n(a)i  )    for  each 

positive   y.   Let   p(u)  =  i   if   m .  <_  n(a)  <  m.  ,  ;  then   p(a) 

tends  to  infinity  as   a  ->-  °°   but 


(4.6)    lim   n(a)  P{C  .  .(a)  >  y  /n  (a) } 

P  (a)     — 

a  -*■   °° 


lim   n(a)  p(a)  P{|D;,(a)|  ^y   /n  ( a )  p  ( a )  1) 

a  ■>  °° 


lim   n(a)  p(a)  P{|D1(a)|  >_  y  /n  (a)  p(a)  2} 

a  ■*•  °° 


2  -1 
<_  lim   n(a)  p(a)      [y    n  (a)  p(cx)  ]    =  0 


We  can  further  choose   p(a)   so  that   p(a)  <  n (a) .    Fix  e  >  0 
By  the  previous  lemma  we  can  choose   y   so  that   y  >  a      and 
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(4.7)  P{|Vk(oO  (CX)  '     >    Y    /F^T}    <    £^    2 


for   all      a.      Let 


A. (a)    =  {max  |v^ (a) |     <    3y    /n(a)     <     |V, (a) | }     ; 
then 


1        V<i   j 


P{    max         Iv.  (a)  |     >_   3y    /n  (a)  } 
i  <  n  (a) 


<    P(  lvn(a)  (a)  I     >    y    ZnTaT} 


n (a) -1 


+         J         P(A     (a)     n    (|V     ,    ,     -    V.  (a)  I     >    2y       /FTT^D     ) 

•  _  i  1  n  \  06  /  1 


The  summation  is  bounded  above  by 


n (a) -p (a) -1 

)      P  {  I  V .  ( a  )  -  V  .    ,  x  (  " )  I  >  Y  ^n(a)  } 
.  L ,         l        l+p ( a)      — 
1=1 


n(a)-p(a)-l 

i|i    P<Ai(a)  H  dVn(a)(a)  -  Vi+p(a)(a)|  >y    /KW>}) 

n  (a)  -1 

+      I  P{|Vn,  ,  (a)  -  V.  (a)  |  >  y   /n(a)  }   . 

/  \        i    \  n (a)        l      — 

i=n(a)-p(a) 

Each  term  in  the  first  and  third  of  these  sums  is  bounded  above 

* 
by   p{C  ,  v  (a)  >  y  /n ( a ) } .   Since   A. (a)   is  measurable  with 

respect  to  the  o-algebra  generated  by   U,  (•<),...,    U.(a)   by  (3.7) 
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p{   max   |v.  (a)  I  >_  3y  /n  ( a )  } 
i  <  n  (a) 


-  P{  ,Vn(a)  (a)  '  -  Y  /n(a)  } 


+  [n(a)  +p(a)]  P{C     (a)  >_  y   /n  (a)  } 


n  (a)  -p  (a)  -1 
+  P(Ai(a)  [P{  |Vn(a)  (a)  -  Vj_   (a)  (a)  |  >  y  /n7^r}+4,(p(a)] 

-2  * 

<_  2ey      +<j)(p(a))  +  [n(a)  +p(a)]  P{C  (a)(a)  >_  y  /n  ( a )  } 


by  (4.7)  and  the  fact  that  the   A. (a)   are  disjoint.   Hence  by 
(4.6)  and  (3.8)  for   a   large 


P{   max   |V.  (a)  |  >  3y  /n(a)  }  <_  3ey 
i  <  n(a)   X 


-2 


PROOF  OF  THEOREM  (4.1).   For   a  =  1:   By  Lemma  (4.5)  the 

sequence   {Y  }   is  tight  on  ID [0,1].   Further   lim      E[Y  (t)]  =  0 

and 


lim   E[Y  (t)2]  =  lim   ( /n (a)  a (a) )  2  E 


a 


ry 


[n(a)t] 


/  m  iou  t  j      \ 
D  (a) 

v   k=l    k    ' 


=  t 


by  (3.5).   Since   Y^t)2  =  [a2  (a)    n(a)]  X  v  [n  ((X)  t  ]  (a)  2  ' 

2 
{Y_(t)}   is  uniformly  integrable  for  each   t   by  Lemma  (4.4). 

By  Lemma  (4.2)   {Y^  (t) ;  0  <  t  <  1}   has  asymptotically  independent 

increments.   The  result  now  follows  from  Theorem  (19.2)  of 

Billingsley  [ 1 ] . 
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5.   EXTENSIONS 

This  paper  presented  a  heavy  traffic  result  for  the 
limiting  distribution  of  the  customer  waiting  time  in  single 
server  queues  with  dependent  interarrival  and  service  times. 
There  are  many  different  ways  in  which  to  model  dependent 
interarrival  and  service  times  in  queues  using  the  EARMA  scheme 
Some  other  ways  are  given  in  Lewis  and  Shedler  [7],   Properties 
of  these  multivariate  exponential  models  will  be  given 
elsewhere . 

Work  is  currently  going  on  in  the  simulation  of  these 
queues.   Some  simulation  results  have  been  reported  for  the 
queue  of  Example  (2.7)  with   3=0.5,   p  =  4   and  various 
values  of   P   and   A   in  Boonsong  [2],   They  seem  to  indicate 
that  the  heavy  traffic  result  approximates  the  average 
customer  waiting  time  fairly  well  for  a  traffic  intensity  of 
.9  and  quite  well  for  a  traffic  intensity  of  .95.   However, 
the  distribution  of   W   in  a  simulation  with   p  =  0.25  and 
traffic  intensity  .995  shows  a  large  amount  of  underdispersion 
relative  to  an  exponential  distribution.   Thus,  it  seems  that 
the  convergence  of  the  waiting  time  distribution  is  very  slow. 

It  is  expected  that  further  simulations  will  give 

more  indication  of  the  validity  of  using  the  heavy  traffic 

result  to  approximate  the  average  customer  waiting  time,  or 

quantiles  of  the  waiting  time  distribution  under  less  than 

heavy  traffic  conditions.   These  results  will  be  reported 

on  elsewhere. 
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